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Ultrahigh energy neutrons and pions are likely to be produced in particle interactions inside 
cosmic ray sources and subsequently decay to neutrinos and other secondary particles [ty — > 



(Vfi),!^ 



e Vpivp^Veiy^. In high magnetic fields of the cosmic acceleration sites, the 



ultrahigh-energy-charged particles may lose energy significantly due to synchrotron radiation before 
decay. We show that for gamma ray bursts in the internal shock model the flux of very high-energy 
antineutrinos (v e ) produced from decaying ultrahigh-energy neutrons can be more than the total 
neutrino flux produced in pion decay depending on the values of their Lorentz factors, luminosities 
and variability times. 

PACS numbers: 98.70.Dk, 98.70.Sa 



O 

(N 

in 



Oh- 
6 

CO 

(N 
> 

O 

in 

(N 

q 

in 

o 
o 



X 



Pierre Auger [U [2| and other extensive air- 
shower arrays [iMl have detected large number 
of ultrahigh-energy cosmic ray events. The ori- 
gin of these cosmic rays are yet to be identified. 
Some of the ultrahigh-energy cosmic ray events 
have been correlated with active galactic nuclei 
(AGN). These events may also come from other 
sources with similar spatial distributions. Gamma 
ray bursts (GRBs) have long been speculated as 
sources of ultrahigh energy cosmic rays [6|, |7|. 
High-energy gamma rays have been detected from 
many GRBs 8] . These gamma rays can be of lep- 
tonic or hadronic origin. More observational data 
may constrain the gamma ray production mecha- 
nisms in the near future. As protons can be ac- 
celerated to 10 21 eV inside GRBs by Fermi mecha- 
nism, they are very attractive candidates for mul- 
timessenger astronomy. The shock-accelerated ul- 
trarelativistic particles lose energy by interactions 
and radiations inside their sources and also during 
propagation. The charged cosmic rays are largely 
deflected by Galactic and interstellar magnetic 
fields and it is difficult to trace back their origin. 
Accelerated protons interact with low energy pho- 
tons and protons inside their sources. Charged and 
neutral pions are generated in this way. Charged 
pions subsequenly decay to high energy secondary 
neutrinos pj, pp — > n X, tt^ 1 — »■ /i ± ^ A ,(i/ M ), fj^ — > 
e ± D fl (iy fl )i' e (D e ). Ultrahigh-energy gamma rays are 
produced from decay of the neutral pions 7r° — >■ 77. 
The high-energy gamma rays travelling from high- 
redshift sources are likely to be absorbed by in- 
frared background and lower-energy photon show- 
ers will be produced in cascade interactions. They 
are found to be useful in probing nearby acceler- 
ation sites [9J . Moreover, ultrahigh-energy pho- 
tons produced by extremely energetic protons in- 
teracting with cosmic microwave background ra- 
diation, during propagation through the interstel- 
lar medium, can act as messengers of the cos- 



mic accelerators [lfj. The ultrahigh-energy sec- 
ondary neutrinos produced in pp and pj interac- 
tions inside the sources are expected to travel to 
us undeflected by Galactic and extragalactic mag- 
netic fields. Neutrino telescopes like IceCube [ll| . 
ANITA [12|] are searching for neutrino signals to 
complement the gamma-ray detectors. Simultane- 
ous observations by various detectors may reveal 
many aspects of the same phenomena. Neutrino 
fluxes from GRBs [TJUl and AGN \fi$Hk have 
been calculated earlier. The neutrino flux from 
Centaurus A has been calculated 2CH22| consider- 
ing P7 and pp processes and found to be detectable 
by IceCube in a few years of operation. If the mag- 
netic field inside the sources is very high then the 
ultrahigh-energy-charged pions and muons lose en- 
ergy significantly before decaying to secondaries. 
As a result their secondary neutrino flux is very 
low. Shock-accelerated protons are expected to 
produce high energy neutrons in various interac- 
tions; these unstable particles subsequently decay 
to leptons and protons. This channel was previ- 
ously considered [23[ to estimate TeV neutrino flux 
from Cygnus OB2 located about 1.7 kpc away from 
us. In this paper we discuss that the neutron decay 
channel of ultrahigh energy neutrino production 
may become more important than the pion and 
muon decay channels in case of cosmic accelerators 
with high internal magnetic fields. In particular, 
we have considered the neutrino production from 
GRBs in the internal shock model through pho- 
topion and neutron decay, including the effect of 
synchrotron energy loss by pions and muons in the 
internal magnetic field of fireball. The values of 
the Lorentz factor, luminosity, and the variability 
time of a GRB determine whether this effect will 
be significant at high energy. In future it would 
be interesting to study this effect in other sources. 
The secondary antineutrinos from neutron decay 
are particularly useful to put upper limit on the 



energy of the cosmic-ray protons originating from 
cosmic acceleration sites of high magnetic fields. 



HIGH ENERGY NEUTRINOS FROM 
GAMMA RAY BURSTS 

We discuss the photopion and neutron decay 
channels of high-energy neutrino production in- 
side GRBs. The high-energy neutrino flux gen- 
erated in pj interactions through pion decay from 
individual GRBs in the internal shock model has 
been calculated earlier 
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In [16j the au- 
thors have derived the neutrino flux from individ- 
ual GRBs using the isotropic photon energy typ- 
ically observed by Swift in the energy range of 1 
KeV to 10 MeV and folded the neutrino flux with 
the distribution functions to incorporate the fluc- 
tuations in the values of the GRB parameters. The 
diffuse neutrino flux from all high and low- lumi- 
nosity GRBs have been calculated in 16| consider- 
ing distribution functions for the GRB parameters. 
GRBs were assumed to follow star formation rate 
and the rate of low luminosity GRBs was taken to 
be 500 times higher than high luminosity GRBs. 
In this paper we have calculated the neutrino flux 
from individual GRBs by expressing it in terms of 
the observed isotropic photon energy as discussed 
in [16] assuming certain values for the luminosity, 
Lorentz factor and variability time of a GRB. We 
have not considered the distribution in the values 
of the GRB parameters which is important in case 
of the diffuse flux from all GRBs. The photon en- 
ergy spectrum from a GRB can be expressed as a 
broken-power law with break at et, r in the source 
rest frame. 
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These photons are interacting with shock ac- 
celerated protons to produce charged and neutral 
pions. The charged pions subsequently decay to 
muons and neutrinos. Muons also decay to neu- 
trinos. The minimum energy of the protons inter- 
acting with photons of energy e^r is 

_1 GeV. 
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In P7 interactions both 7r° and tt + can be produced 
with equal probabilities. 7r + gets, on average 20% 
of the proton energy and if the final-state leptons 
share the pion energy equally then each neutrino 
carries 5% of the initial proton energy. The first 
break in energy in the neutrino spectrum, E v \, is 
due to the break in the photon spectrum at Cb r - 
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The fire ball Lorentz factor Taoo = T/300, photon 
luminosity L. 



7,51 



L 7 /(10 51 ergs /sec), variability 
time £t,,-3 = t v /(10~ 3 )sec are the important pa- 
rameters in our calculations. The total energy to 
be emitted by neutrinos of energy E v can be ex- 
pressed as, 
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where E l *° is the total isotropic energy of the emit- 
ted gamma-ray photons in the energy range of 1 
KeV to 10 MeV; it is available from observations. 
ee,e_B ~ 0.3 are the energy fractions carried by 
electrons and the magnetic field respectively. /„. 
is the fraction of the energy of a proton going 
to pion production inside a GRB fireball. The 
relativistic electrons produce the photons by syn- 
chrotron radiation, so four orders of magnitude in 
photon energy corresponds to two orders of magni- 
tude in the energy of the radiating charged leptons. 
We have assumed that the electrons are radiating 
all their kinetic energy as observed photon energy 
E™°. The normalization constant of the power law 
electron flux is proportional to E lso / k. The differ- 
ential spectral index of the electron spectrum is 
assumed to be —2.5 which corresponds to K = 1.8. 
The total energy carried by the differential flux of 
neutrinos has been approximately related to the 
total energy carried by the differential flux of elec- 
trons in Eq.(Q|. The first break [Eq.Q] in the neu- 
trino spectrum due to the break energy in the pho- 
ton spectrum is contained in the expression of / w 
[16J . The charged pions and muons lose energy due 
to synchrotron radiations in high magnetic fields 
inside the sources before decaying to neutrinos. As 
a result, the neutrino flux is depleted at very high 
energy. The neutrino flux from photopion decay 
in the source rest frame is 
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(5) 
Pion-cooling energy is 10 times higher than muon- 
cooling energy. We would be overestimating the 
total neutrino flux if we use the pion-cooling en- 
ergy to derive the second break energy in the neu- 
trino spectrum. The muon-cooling break E* can 
be expressed as a function of GRB parameters. 
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The fraction of fireball proton energy lost 
neutron production is, 



f n (E p ) = 0.8/q 
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are discussed in [16| . The average value of the frac- 
tion of energy transferred from the proton to the 
neutron is assumed to be (x p ^. n ) = 0.8. The neu- 
tron spectrum can be expressed in terms of E. 
and the equipartition parameters e e , cb- 
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2 corresponds to the assumption that charged and 
neutral pions are produced with equal probabili- 
ties. The GRB is assumed to be at a redshift of 
0.8 in all our calculations. Duration of the burst 
is assumed to be 5 sec. The break in the pho- 
ton spectrum is assumed to be at e{, r = O.SMeV 
in the source rest frame. We have assumed the 
spectral indices of the photon flux 71 =1.5 and 
72 = 2.2. The ultrahigh energy neutrons decay 
(n — > p + e~ + V e ) to v e . Their rest frame life- 
time is t„ = 886 seconds and decay mean free 
path 10(E n /EeV) kpc. The maximum energy of 
the neutrons producing v e depends on the distance 
D s of the source. If the source is at a large dis- 
tance then all the neutrons are expected to decay 
during propagation. The v e flux expected on earth 
at energy E p from decaying neutrons ( mass m n ) 
of flux (dN n /dE n ) [23J is 
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where eo is the mean energy of v e in the neutron 
rest frame. The maximum energy of the shock ac- 
celerated protons E p _ max in the GRB fireball can 
be 10 12 GeV. The maximum energy of the sec- 
ondary neutrons is E n ^ max — 0.8E Ptmax , assum- 
ing 80% of the proton's energy goes to a neu- 
tron in P7 interactions. The maximum energy 
of the antineutrinos is Q/m n = 0.001 times the 
maximum energy of the decaying neutrons, where 
Q = m n — m p — m e . The observed antineu- 
trino/neutrino flux on earth is 



3. respectively. Figure 2. shows for high values of 
luminosity this effect is more significant. The num- 
ber of neutrino events expected from a GRB in Ice- 
Cube is very low. The neutron decay channel can 
give only 10 -5 muon events in Km 2 area above 1 
TeV energy from a GRB at redshift 0.8, with lumi- 
nosity L 1 — 10 53 erg/sec, Lorentz factor T — 300, 
duration 5sec, and variablity time t v — 10ms. The 
GRB has been assumed to occur at a zenith angle 
of 120° with respect to an observer at the south 
pole. The number of total neutrino events con- 
tained inside 2fcm 3 volume of ice is expected to 
be 10~ 5 from the neutron decay channel and 10~ 4 
from the pion decay channel. Above 10 16 or 10 17 
eV we expect the neutron decay channel to dom- 
inate over the pion decay channel depending on 
the values of the GRB parameters. Neutrino tele- 
scopes like ANITA may be useful for detecting the 
neutrino fluxes in this energy range. 



CONCLUSION 

Multimessenger approach to explore the cosmic- 
ray universe has gained much importance due 
to the successful operations of the cosmic-ray, 
gamma-ray and neutrino detectors. In this pa- 
per we have discussed that the antineutrino flux 
produced in decay of ultrahigh-energy neutrons 
may become more important than the neutrino 
flux from photopion decay depending on the lu- 
minosity , Lorentz boost factor, varibility time of 
a GRB at very high energy. In future one may 
consider other cosmic acceleration sites with high 
magnetic fields to explore the importance of the 
neutron decay channel at very high energy. 
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Figures 1-3 show our calculated neutrino fluxes 
from photopion decay (dashed line) and neutron 
decay (solid line) in the observer's frame. In Fig.l. 
the Lorentz factor has been varied from 100 to 
500 assuming luminosity 10 51 erg/sec and variabil- 
ity time 10ms. For small values of Lorentz factor 
and variability time the neutron decay channel can 
become more important than the photo-pion decay 
channel at very high energy, as shown in Figs.l and 
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FIG. 1: Neutrino flux from pion decay (dashed line) and neutron decay (solid line) from GRB with T — 100, 200, 
300, 500 (top to bottom), L = 10 51 erg/sec, t v = 10 ms. 










FIG. 2: Neutrino flux pion decay (dashed line) and neutron decay (solid line) from GRB with L = 10 52 , 10 51 , 
10 50 erg/sec (top to bottom), F = 300, t v - 10 ms. 
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FIG. 3: Neutrino flux from pion decay (dashed line) and neutron decay (solid line) from GRB with t v — 1, 5, 
20 ms (top to bottom), F = 300, L= 10 51 erg/sec. 
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